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Abstract 

The phenomenology of a system of two coupled quadratic maps 
is studied both analytically and numerically. Conditions for synchro- 
nization are given and the bifurcations of periodic orbits from this 
regime are identified. In addition, we show that an arbitrarily large 
number of distinct stable periodic orbits may be obtained when the 
maps parameter is at the Feigenbaum period-doubling accumulation 
point. An estimate is given for the coupling strength needed to obtain 
any given number of stable orbits. 



1 Introduction 

A multistable system is one that possesses a large number of coexisting at- 
tractors for a fixed set of parameters. There is ample evidence for such 
phenomena in the natural sciences, with examples coming from neurosciences 
and neural dynamics 0-0, optics || 0, chemistry [TO | [[31J [ l2| , condensed 



matter []n| and geophysics [HI]] . Multistability also seems to be an essential 



complexity-generating mechanism in a large class of agent-based models [15 

In view of this, it is important to identify the dynamical mechanisms 
leading to multistability and, in particular, to construct simple models where 
this phenomenon might be under control. The first mathematical result in 
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this direction was obtained by Newhouse ]If| ]I7J ]Tj| who proved that, 
near a homoclinic tangency, a class of diffeomorphisms in a two-dimensional 
manifold has infinitely many attracting periodic orbits (sinks), a result that 
was later extended to higher dimensions It has also been proved [20 



that, in addition to infinitely many sinks, infinitely many strange attractors 
exist near the homoclinic tangencies. The stability of the phenomena under 
small random perturbations has been studied |2l]] |22|| . 

A second dynamical mechanism leading to multistability is the addition 
of small dissipative perturbations to conservative systems. Conservative sys- 
tems have a large number of coexisting invariant sets, namely periodic orbits, 
invariant tori and cantori. By adding a small amount of dissipation to a con- 
servative system one finds that some of the invariant sets become attractors. 
Not all invariant sets of the conservative system will survive when the dissi- 
pation is added. However, for sufficiently small dissipation, many attractors 
(mainly periodic orbits) have been observed in typical systems [|23| |24[ |25| 



The problem of migration between attractors and their stability in multiple- 
attractor systems has also been studied by other authors [^6| P7| . Most of 
results are based on numerical evidence. However, using the techniques of 
deformation stability [^HJ [^9| J3(J [|3T| some rigorous mathematical results 
p2| may be obtained. 

Finally, it has been found recently that, for parameter values near the 
Feigenbaum period-doubling accumulation point, quadratic maps coupled by 
convex coupling may have a large number of stable periodic orbits. This is 
one of the phenomena we study in detail in this paper. The emphasis on 
quadratic maps near the Feigenbaum accumulation point has a motivation 
close to the idea of control of chaos [0 P5|| . The typical situation in control 
of chaos, is that of a strange attractor with an infinite number of embedded 
periodic orbits, all of them unstable. These orbits are then stabilized by 
several methods. If, instead of a large number of unstable periodic orbits, 
one has, for example, a large number of sinks, the controlling situation would 
seem more promising and robust, because the control need not be so accurate. 
It would suffice to keep the system inside the desired basin of attraction. At 
the period-doubling accumulation point the Feigenbaum attractor, because of 
the properties of the flip bifurcations, coexists with an infinite set of unstable 
periodic orbits. By coupling, as we will show, an arbitrarily large number of 
orbits may become stable. 

The existence of a large number of stable periodic orbits for just two 
coupled quadratic maps, provides a simple model where multistability is well 
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under control, in the sense that not only the nature of the phenomenon is 
completely understood as one may also compute the range of parameters that 
provides any desired number of stable orbits. This should be contrasted, for 
example, with concrete models for the Newhouse phenomenon [[36] . 



Rather than merely focusing on multistability, we also study the phe- 
nomenology of two coupled quadratic maps, in particular the bifurcations 
of periodic orbits and the regime of synchronization.. The stabilization of 
orbits in the coupled system is similar to that obtained in higher dimensional 
coupled map lattices ]37j with the exception that, due to the restricted di- 
mension of the phase space, the types of bifurcations are different in our 
system. The results concerning the multistability phenomenon at /i = /i 



oo 



also considerably extend, and also correct, some imprecise statements in |33 



2 Coupled quadratic maps 

Coupled map lattices (CML) are discrete dynamical systems generated by 
the composition of a local nonlinearity and a coupling. The phase space of 
the CML considered in this letter is the square [— 1, l] 2 and the dynamics is 
generated by the map F e defined as follows. Given a point (x, y) G M., its 
image by F e , denoted (x, y) is given by 

x=(l-e)f(x) + ef(y) m 
y = (l-e)f(y) + ef(x) {L> 

where < e < | , f(x) = 1 — fix 2 and < ji < 2. 

The map / maps [—1,1] into itself. Therefore, the convex combination 
in (|^) ensures that F e ([—1, l] 2 ) C [—1, l] 2 and the dynamics is well-defined. 
We denote the orbit issued from the initial condition (x, y) by the sequence 
{(x t ,y t )} t ^n, that is to say, (x°,y°) = (x,y) and (x t+1 ,y t+1 ) = F e (x t ,y t ) for 
all t E N. 

For the sake of simplicity, we will often employ the variables s = and 
d = ^r~. The previous notation of orbits also applies to these variables for 
which relation ([!]) becomes 

' s_= l-fi 2 (s 2 + d 2 ) (2) 
d = —asd 

where a = 2yu(l — 2e). 

Finally, note that the dynamics commutes with the symmetry a(s,d) = 
(s, —d) or a(x,y) = (y,x) in the original variables. 
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3 Synchronization 



If d l = 0, then d t+1 = and s t+1 = /(s*). In this case, the orbit is said to be 
synchronized (from t on). More generally, an orbit is said to synchronize if 
linit^oo rf* = and if all orbits synchronize, then we say to have synchroniza- 
tion of the map. Synchronization is the simplest dynamical regime exhibited 
by two-dimensional CML. 

To determine a sufficient condition for synchronization in our system, we 
note that for any orbit, one has |s*| < 1 for all t G N. It follows from (0) 
that the condition \a\ < 1 ensures an exponential decay of and hence 
synchronization.. Since a > 0, the condition \a\ < 1 is equivalent to the 
following ones (see Figure 1). 

t = \ ° r e< \ and " < ^ l(e) := 2(T^) ' 

From now on, we assume that e < |. The condition /i < /xi(e) is not 
necessary. Indeed, if for instance [i < |, then / has an attracting fixed 
point in [—1,1], and one can prove that synchronization occurs. When e is 
sufficiently small, this happens even though \i > /xi(e). 



4 Non-synchronized period-2 orbits 

Starting from synchronization and modifying the parameters, non-synchronized 
(periodic) orbits appear from bifurcations of synchronized (periodic) ones. To 
understand this phenomenon, as well as the bifurcations of subsequent orbits, 
we now study analytically the periodic orbits of period 1 and 2. Let (s*,0) 
be the synchronized fixed point denoted in (s, <i)-variables. It exists for any 
values of the parameters and F € has no other fixed point in [—1, l] 2 . In (s, d)- 
variables, the Jacobian of F € at this fixed point is diagonal. One eigenvalue 
is f'{s*) and the corresponding eigendirection is the diagonal d = 0. The 
other eigenvalue is (1 — 2e)/'(s*) and the corresponding direction, orthogonal 
to the diagonal is referred as the anti-diagonal. 

The condition f'(s*) = —1 then determines a period-doubling bifurcation, 
which is of co-dimension 1 if e > 0. This is the well-known period-doubling 
bifurcation of / which creates a synchronized period-2 orbit of F t . 

Moreover, one checks that the derivative f'(s*) is negative for any //. 
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Hence the conditions e > and (1 — 2e)f'(s*) = — 1, i.e. 

3 -4e 



e > and /i = /12(e) 



4(1 -2e) 2 ' 

determine another co-dimension 1 period-doubling bifurcation of the syn- 
chronized fixed point. Indeed the conditions of the corresponding bifurca- 



tion theorem (see e.g. [PS[| ) are satisfied when the curve /1 = /12(e) is crossed 
upward. 

The period-2 orbit created at this bifurcation is non- synchronized and 
symmetric. To show this, denote by (s\,di) and (32,^2) its components. 
Since the multiplier (1 — 2e)f'(s*) is negative and the bifurcating direction 
is the anti-diagonal, we have did?, < (sufficiently close to the bifurcation). 
Because of the a symmetry, the map F e also has a period-2 orbit with com- 
ponents (si, — di) and (s 2 ,— c^)- Consequently if a(s%,di) 7^ (£2,^2), the 
system would have two periodic orbits created by a co-dimension 1 bifurca- 
tion. This is impossible by the unicity in the bifurcation theorem. Therefore, 
sufficiently close to the bifurcation, we have a(si,di) = (52,^2) and di / 
which is the desired conclusion. 

By continuity in the parameters of F e , sufficiently close to the bifur- 
cation, this symmetric orbit is stable with respect to perturbations in one 
direction (the anti-diagonal direction at the bifurcation) and since f'{s*) < 
(1 — 2e)f'(s*) < —1, it is unstable in the direction orthogonal to the latter. 

The bifurcations will now be computed. The orbit with a(si,di) = 
(s2,d 2 ) and di ^ exists for any /1 > /i 2 (e) and is the unique (up to time 
translations) period-2 non-synchronized symmetric orbit of F 6 in [— 1, l] 2 . 

Computing the corresponding Jacobian, one obtains the equation for the 
multipliers 

A 2 - 2 [(1 - 2e) (S 2 - D 2 ) + 2e 2 S 2 ] A + (1 - 2ef (S 2 - D 2 ) 2 = 

where 




l-2e 



_ A /4 At (l-2 £ ) 2 -(3-4 £ ) 
l-2e 

Direct calculations show that, if e > 0, the multipliers, say A_ and A + , 
have zero imaginary part iff /i < ^z^s, ■ Under this condition, we have 
< A_ < 1 if fj, > /i 2 (e) and A+ > 1 iff /x 2 (e) < fi < /13(e) := 4(1 ^ 2e)2 ■ (The 
inequality /12(e) < /13(e) indeed holds if < /1 < 2 and < e < |, see Figure 
1.) 
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Consequently, by increasing /x, the symmetric orbit suffers an inverse 
pitchfork bifurcation at \i = [13(e)- This bifurcation is generic for a symmetric 
orbit in a system with symmetry |38| and the conditions of the bifurcation 
theorem hold when the curve \i = /43(e) is crossed upward. 

This bifurcation creates two non-symmetric period-2 orbits (one orbit and 
its symmetric). We have checked that these orbits exist for any e > and 
/i > /43(e). For e = 0, their components are combinations of a fixed point of 
/ and the components of a period-2 orbit. 

When the imaginary part of A_ and A + is not zero, we have |A_| = |A+| < 
1 iff /j, < //4(e) := 4(fz§t)2- (Once again, if < /4 < 2 and < e < |, the 
inequality //3(e) < /44(e) is satisfied, see Figure 1.) 

The symmetric orbit is thus stable in the interval /42(e) < /t < //4(e). If 
e > and the curve /t = /14(e) is crossed upward, this orbit suffers a Hopf 
bifurcation creating a locally stable invariant circle. A numerical calculation 
shows that the latter is destroyed when \x is sufficiently large or when e is 
sufficiently small. Obviously, if e = 0, it does not exist and the bifurcation at 
/t = /44(C)) = | which is a period-doubling bifurcation of / creates a period-4 
orbit. 

Note that invariant circles in two-dimensional CML resulting from the 
destabilization of a symmetric orbit and their normal form had already been 
reported in |39|]. In that work, the system is also defined by ([[]), but the 
local map is fix) = ax(l — x) and e may be larger than |. Figure 2 shows 
an example of the phenomenology described above. Numerically, it is more 
convenient to follow the orbits from e = to increasing values of the coupling. 
In this picture, as well as in the following ones the map parameter is fi = /too 
(the accumulation point of the period-doubling cascade). In Figure 2, from 
e = (the circle) to e = 0.039 (the point labelled 1), the symmetric orbit is 
unstable. The figure also shows the invariant circle for e = 0.03. Between 
the points 1 and 2, the symmetric orbit is stable. At point 2, the pitchfork 
occurs, the symmetric orbit becomes unstable and the non-symmetric orbits 
are created. Finally, the point 3 (e = 0.182) corresponds to the collapse on 
the synchronized fixed point. 

5 The phase opposition orbits 

The previous phenomenology is not restricted to small periods but extend to 
any power of 2. In particular, the synchronized period-2 p orbit may destabi- 
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lize to create a symmetric (non-synchronized) orbit of twice the period. 

Given p e N, let {sj}i<j<2P be the components of the period- 2 P orbit of /. 
The points {(sj, 0)}i<j<2P are the components of the synchronized period- 2 P 
orbit of F e . By the chain rule and since each Jacobian at (s*, 0) is diagonal, 
the corresponding multiplier along the anti-diagonal direction is 

2 p 

(1-26)^11 

i=i 

The condition that this multiplier equals —1 determines, if e > 0, a co- 
dimension 1 period-doubling bifurcation. Applying the reasoning of the pre- 
vious section to each component (sj,0), we conclude that this bifurcation 
creates an orbit with the property <r(s i+2P , d t+2P ) = (s*, d l ) and d l ^ for all 
t e N, which is called a phase-opposition period- 2 P+1 orbit. 

Since (1 — 2e) 2P < 1, this bifurcation occurs only if the bifurcation along 
the diagonal direction has occurred (the local period-doubling bifurcation of 
{si}i<i<2p)- In other words, the phase opposition period-2 p orbit exists only 
if the synchronized period- 2 P orbit does. Moreover it follows from Figure 3 
that, at least for ji = /i^, the phase opposition period- 2 P+1 orbit exists only 
if the phase opposition period- 2 P orbit does. This is confirmed analytically 
for the period— 4 orbit whose existence condition is the instability of the 
synchronized period— 2 orbit in the anti-diagonal direction. One obtains 

1 

/!>/! ( 4)(e) = l + 4(i _ 2g)2 , 

and /i2(e) < A*(4)(e) if < e < | and /^(e) > 2 and //(4)(e) >2if|<e<| 
(see Figure 1). 

Furthermore, a numerical calculation at fi = /x^ , reported in Figure 3, 
shows that the succession of bifurcations of a phase opposition orbit does 
not depend on the period. On this picture, we have plotted the values of e 
for the Hopf bifurcation, the pitchfork bifurcation and the period-doubling 
bifurcation creating the orbit, versus the power of the period. For each 
period, the phenomenology is identical to that described in the previous 
section, with an adequate change of scale in e. In addition, the picture shows 
that several phase opposition orbits may be stable for e > fixed. This 
stabilization is an effect of the coupling that will be discussed below. 

Finally, since the phase opposition orbits are the first orbits to appear 
when the parameters are varied from synchronization and since the first such 
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orbit that is created is of period 2, it follows that a necessary and sufficient 
condition for synchronization is fj, < //2(e), the condition for the existence of 
the latter. 

6 The non-symmetric orbits 

We now analyze the existence and the stability of other period-2 p orbits for 
A* = Aioo- Our interest for this value of \i is that the scaling properties of 
/ are reflected on scaling laws for the periods and values of e at which the 
bifurcations occur (see Figure 3 and 8). We only consider the orbits which for 
e = have the same period on projection to both axis x and y. These orbits 
are followed numerically when e increases and are referred using the phase 
shift of their components at e = 0. For fj, = /i^, the map / has a period- 2 P 
orbit for each pGN, whose components for p up to 5 are shown in Figure 5. 
In this picture, the numbers reflect the order in which the components are 
visited and the tree structure represents the origin of each component in the 
bifurcation cascade. An important notion is the dyadic distance 5 between 
the components of an orbit. 5 is the number of steps one has to go back in 
the bifurcation tree to meet a common component. 

The dyadic distance is used to characterize the families of periodic orbits 
that we are considering. For instance, the coordinates of each component of 
a synchronized orbit are at distance 0, those of a phase-opposition orbit are 
at distance 1. Accordingly, when we speak of distance— k orbit we refer to 
the dyadic distance of the coordinates of its components. 

For any 5 > 1, there are 2 5 ~ l different orbits with distance 5 which have 
coordinates out of phase by 2 p - 5 + a 1 2 p ~ m + a 2 2 p ~ s+2 + ■ ■ ■ + a d _ 1 2 p ' 1 steps, 
with ccj e {0, 1}. The distance of a period- 2 P orbit is at most p (p e N). 
For e = 0, the only symmetric orbits are those at distance and I. This 
property is preserved for e > as shows Figure 6 for 5 = 2. The succession 
of bifurcations of orbits with distance 5 > 2 should then differ from those with 
distance 1. The differences are seen in Figure 6 which shows the evolution of 
the eigenvalues. For e = 0, the orbit is unstable. When e increases, it suffers 
two collisions with orbits of twice the period when the eigenvalues cross —1 
and then becomes stable. (When decreasing e, these collisions would be 
period-doubling bifurcations.) If e increases further, the orbit collides with 
an unstable one of the same period in a saddle-node bifurcation when the 
larger eigenvalue reaches 1. For larger values of e, the orbit does not exist. 
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The unstable orbit with which it collides is the one that at e = has period 
2 P in one projection and 2 P_1 in the other. For higher dyadic distances, the 
overall variation of the eigenvalues is similar to the 5 = 2 case. Figure 7 
shows a typical example of these phenomena for the case 6 = 2. Between 
e = and the point labelled 1 in the figure, the orbit is unstable. The point 
1 corresponds to the smaller eigenvalue crossing -1 (see Figure 6). Therefore, 
between the point 1 and 2, the orbit is stable. It disappear at the point 2 
when it collides with an unstable orbit of the same period. 



7 Multistability 

We have seen that the coupling stabilizes the orbits with distance larger 
than at fi = /ioo. There are indeed two mechanisms responsible for this 
stabilization.. The determinant of the Jacobian of a period-2 p orbit is 

(i- 2e ) 2p (-2 / . 0O ) 2p+i n^- 

i=i 

The term (1 — 2e) 2P coming from the coupling decreases when p increases. 
However, there is yet a second stabilizing mechanism. Denote by T(e) the 
remaining factor in the determinant 

T(e) = (-2 / u 00 ) 2P+1 JJxi^/i. 

i=i 

Without coupling, T(0) is simply the square of the multiplier of / for the 
periodic orbit. From the properties of the Feigenbaum - Cvitanovic functional 



equation it follows [33] that this factor converges to a fixed value around —1.6 
when p increases. The coupling however, changes the position of the orbit 
components in such a way that this factor also decreases. It is the combined 
action of this decrease with the contraction of the coupling that brings the 
eigenvalues into the interior of the unit circle and stabilizes the orbits. 

For small e there is a simple geometrical interpretation for the variation 
of T(e). The reason why in the one dimensional map the product Y\a=i x i 
remains constant, when p grows, is because each time the period doubles, 
the doubling in the number of factors greater than one is compensated by 
the fact that the component of the orbit closest to zero approaches zero a 
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little more. For the unstable orbits along the period-doubling chain, the 
orbit components closest to zero alternate on each side of the origin. The 
contracting effect of the convex coupling tends to bring the orbits back in the 
period-doubling hierarchy. Therefore, because the component closest to zero 
has to move across the origin for the orbit to approach the one with half the 
period, this implies that the product of the coordinates is going to decrease. 
The greater the dyadic distance between the orbit projections on the axis, 
the greater will be the perturbation that the original (one-dimensional) orbits 
suffer. Therefore one expects the contracting effect in T(e) to increase with 
the dyadic distance. This effect is quite apparent on Figure 8 which shows the 
stabilizing and destabilizing lines for orbits with distance from 1 to 4. The 
shift downwards of the stable regions for successively larger dyadic distances 
implies that the smaller e is, the larger the number of distinct stable orbits 
that are obtained. An accurate numerical estimate of the number of distinct 
orbits is obtained by computing the derivative D (p, 5) — d (j^ =1 x iVi) I ® e 
at e = for each p and dyadic distance 5. Actually this derivative provides 
an accurate estimate of T (e) itself, because this one varies almost linearly 
with e for most of the stable range of the orbits. On Figure 9, the scaling 
properties, when p grows, of this derivative are shown. From these results 
one computes 

log(-D(p,5))=g(5)-2v +1 

with 

g (5) ~ 0.9075 + 3.475 

Notice that in Figure 9 there is more than one data point for each pair (p, 5) 
which correspond to non-equivalent orbits with the same dyadic distance. 
Two other useful estimates are: 

• the value of the smallest e parameter that stabilizes an orbit of dyadic 
distance 5 equal to the power p 

log eg = -!-535p- 1.547 

• the value of the largest e parameter for which a 5 = 1 orbit is stable 

logcg^ = -0.679p - 1.235 
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From this, one obtains the result that at least k distinct stable orbits are 
obtained if 

< e < exp (-0.99 - 1.22k) 

k is only a lower bound on the number of distinct stable periodic orbits, be- 
cause here we have studied only orbits with the same period under projection 
in the two axis. 

In conclusion: for sufficiently small e an arbitrarily large number of dis- 
tinct stable periodic orbits is obtained. However, for any fixed e, it is an 
arbitrarily large number that is obtained, not an infinite number. Most or- 
bits either synchronize (and are then unstable) or disappear as e grows. As 
a result, a reasoning based on the implicit function theorem, as used in 



is misleading. Given a sequence of orbits of different periods, even if they 
remain as orbits for a small perturbation, that does not mean that their 
(smallest) periods remain distinct. 
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Figure 1: Synchronization curve \x = /J>i(e) and other bifurcations curves, 
/^(e) is the period-doubling of the synchronized fixed point (i.e. birth of the 
phase opposition period-2 orbit), /^(e) is the pitchfork of the phase opposition 
period-2 orbit and //4(e) the Hopf bifurcation of the same orbit. /i(4)(e) is the 
period-doubling of the synchronized period-2 orbit (i.e. birth of the phase 
opposition period-4 orbit) 
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Figure 3: Bifurcation values of e for the phase opposition orbits. From top 
to bottom, birth, pitchfork and Hopf bifurcation 
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Figure 6: e— evolution of distance— 2 orbits and their multipliers at ji — ji, 
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Figure 8: Stability lines of orbits from distance 1 (top) to 4 (bottom) 
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Figure 9: Scaling properties of D (p, 5) 
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